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The structural properties of square-shoulder fluids are derived from the use of the rational 
function approximation method. The computation of both the radial distribution function 
and the static structure factor involves mostly analytical steps, requiring only the numerical 
solution of a single transcendental equation. The comparison with available simulation data 
and with numerical solutions of the Percus— Yevick and hypernetted-chain integral equations 
shows that the present approximation represents an improvement over the Percus- Yevick 
theory for this system and a reasonable compromise between accuracy and simplicity. 
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1. Introduction 



It is well known that the inclusion of attractive interactions in the intermolecular 
potential used to describe a fluid is crucial to obtain a liquid-vapour transition. 
Perhaps the simplest model accounting for this fact is the square-well (SW) fluid 
in which the interaction potential is given by 

{co, r < a, 
-e, a < r < Act, (1) 
0, r > Act, 



where r is the distance, a is the diameter of the hard core, e > is the well depth 
and (A — l)<r is the well width. The thermodynamic properties of the SW fluid 
only depend on three dimensionless parameters, namely the packing fraction r] = 
(7r/6)po" 3 {p being the number density), the reduced temperature T* = ksT/e (ks 
and T being the Boltzmann constant and the absolute temperature, respectively) 
and the width parameter A. Due to the combined assets of relative simplicity and 
'realistic' features, the SW fluid has been studied thoroughly using both theoretical 
approaches and simulations (see, for example, [l-20] an d references therein). 
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Another closely related interaction potential is the 'square-shoulder' (SS) poten- 
tial, that has also been considered in the literature. It reads 



oo, r < a, 

<hs(r) = { e, o < r < Act, (2) 
0, r > Act. 



This purely repulsive potential, apparently considered first by Hemmer and Stell 40 



years ago [2lJ, has been the subject of recent attention 22H4ll]. On the one hand 



it may lead to an isostructural solid-solid transition [22( , to a fluid-solid coexis ting 



line with a maximum melting temperature [231 ] . to unusual phase behaviour [2414261] 
and to a rich variety of (self-organised) ordered structures [27-30]. On the other 
hand, it is the simplest model belonging to the class of core-softened potential 
models for fluids that have been used to study a number of interesting substances 



solutions of electrolytes [34| . 



such as water [3JJ], metallic systems [32], colloidal suspensions [33J and aqueous 



It should also be noted that, as in the SW case, the thermodynamic properties of 
the SS system only depend on the packing fraction 77, the reduced temperature T* 
and the width parameter A. Further, the SS potential becomes equivalent to a hard- 
sphere (HS) interaction of diameter a in the high-temperature limit (T* — > 00) or 
in the limit of vanishing shoulder width (A — > 1), and to an HS interaction of 
diameter la in the low-temperature limit (T* — > 0). 

Despite the simplicity of the SS potential, so far no exact results for the ther- 
modynamic or structural properties of this system are available. Further, even in 
the simplest approximation within the integral equation approach for the study 
of liquids, namely the Percus-Yevick (PY) closure for the Ornstein-Zernike (OZ) 
integral equation, no analytical results have been derived for the SS fluid, Very 
recently, Zhou and Solana [3^] have reported simulations and theoretical results 
based on a bridge function approximation to close the OZ equation, while Guillen- 



Escamilla et al. 40[] have also presented simulation results and a parametrization 
of the direct correlation function which quantitatively agrees with the numerical 
solution of the OZ equation within the PY approximation. 

Several years ago two of us Q derived approximate analytical results for the 
structural properties of the SW fluid using a methodology that has proved useful 
for many other systems [12] • Exploiting the fact that this methodology does not 
make explicit use of the positive character of the well depth, the major aim of this 
paper is to extend it to the SS system and compare the resulting structure with 
the simulation data available to our knowledge in the literature, namely those just 
mentioned by Zhou and Solana [3^ ] and Guillen-Escamilla et al. j^, as well as 
the data in the paper by Lang et al. 0]. Also, given the fact that our approach 
represents an alternative analytical route to the integral equation approach for the 
structural properties of fluids, we will further assess its merits for the SS fluid by 
comparing our results with those we have obtained from the numerical solution of 
the OZ equation both with the PY and the hypernetted-chain (HNC) closures. 

The paper is organised as follows. In order to make it self-contained, in the next 
section we present the derivation of the structural properties of the SS fluid using 
the methodology that we refer to as the Rational Function Approximation (RFA). 
This is followed in Section H by a comparison of our analytical approximation 
and of the numerical solution of the PY and HNC equations for the SS potential 
with the available simulation results. The paper is closed in Section [4] with further 
discussion and concluding remarks. 
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2. Radial distribution function of the square-shoulder fluid 

In this section, we present our proposal for the structural properties of the SS fluid. 
It follows very closely the parallel derivation for the SW fluid as presented in [42j ]. 

2.1. Physical requirements 



As in previous works [43J], it is convenient to consider the Laplace transform of 
rg{r), where g{r) is the radial distribution function (rdf), namely 



oo 



G{s) = / dre- sr rg{r) (3) 



o 

and the auxiliary function ^(s) defined through 

G(s) = S —— T . (4) 

The choice of G(s) as the Laplace transform of rg(r) and the definition of ^(s) 
from Equation @ are suggested by the exact form of g{r) to first order in density 

i- 

Since g(r) = for r < a while g{& + ) = finite, one has 

g (r) = 9(r - a) [g{a + ) + g'(a + )(r - a) + ■ ■ ■ ] , (5) 

where g'(r) = dg{r)/dr and (x) is the Heaviside step function. The foregoing 
property imposes a constraint on the large-s behaviour of G(s), namely 

e° s sG(s) = ag(a + ) + [g(a+) + ag'(a + )] s" 1 + 0(s~ 2 ). (6) 

Therefore, lim< j _ s . 0O e SCT sG(s) = ag(a + ) = finite or, equivalently, 

lim s~ 2 ^(s) = -r— ^ = finite. (7) 

On the other hand, according to the definition of the (reduced) isothermal com- 
pressibility 

dp\ 

dp) 



X = kBT 

POO 

= 1 + 247/cr" 3 / dr r 2 [g(r) - 1] , (8) 
J o 



T 



it follows that 



X = 1 - 247]a 3 lim ^~ 



/ dre~ sr r \g{r) - 1] 
Jo 



= 1 - 247?ct- 3 lim — \G(s) - s~ 2 ] . (9) 
s-s-o ds 

Since the isothermal compressibility x is a ls° finite, one has drr 2 [g(r) — 1] = 
finite, so that the weaker condition J °° drr [g(r) — 1] = lim s _ s .o[G'(s) — s~ 2 ] = finite 
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must hold. This in turn implies 
vp( s ) = -p + pas- \pa 2 s 2 + ( ±p<T 3 + s 3 - (^Ipcr 3 + i-^ as 4 + 0(s 5 ). (10) 
Equation (|4j) can be formally rewritten as 

oo 

= -i E ^ [-*(*)r n (ii) 

n=l 

Thus, the rdf is given by 

1 °° 

S (0 = ^ Vn (r - na) 6 (r - na) , (12) 

71=1 

with 

^ n (r) = -£- 1 {s[-*(s)]- n }, (13) 

/Z" 1 denoting the inverse Laplace transform. 

So far, Equations ©, JTU]) and (fT2j) apply exactly to any interaction potential 
having a hard core at r = a. This includes, among other models, the HS, SW and 
SS potentials. 

For the SS potential given in Equation ([2]), G(s) must reflect the fact that g(r) 
is discontinuous at r = Xa as a consequence of the discontinuity of the potential 
4>ss( r ) an d the continuity of the cavity function y{r) = exp [4>ss( r ) /^bT] g(r). This 
implies that G(s), and hence ^f(s), must contain the exponential term e^( A ~ 1 ) crs . 
This manifests itself in the low-density limit, where the condition limp^o y{ r ) = 1 
yields 



lim g(r) 
p->0 



'0, r < cr, 
e~ l / T * , a < r < Xa, (14) 
1, r > Act, 



lim G(s) = e-^ V -I±^ + fl - e~ 1/r *) e~^i±^. (15) 
Therefore, from Equation (HJ we have 

limtf(a) = — [ e - 1/T *(l + s) + e- (A - 1)s (l-e- 1/T *)(l + As)l~ 1 , (16) 

p->0 Z7T L 

where we have taken, without loss of generality, a = 1. This means that in what 
follows all distances are measured in units of the hard-core diameter a. 

Equations (|ZJ), (fit)]) and (fT6j) are basic exact properties that any reasonable ap- 
proximation for ^(s) must fulfill. 



2.2. The rational function approximation 



In the spirit of the RFA method that has been used in other cases |42j , and following 
the same rationale as for the SW fluid 0], the simplest form that complies with 
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Equation §J§ and is consistent with Equation (116j) is 

*00 



1 — 12?7 + + A 2 s 2 + ^ 3 s 3 
2^ 1 - £ + B lS + e-( A " 1 ) s (S + B 2 s) 



(17) 



where the coefficients Bq, B±, B 2 , A\, A 2 and ^3 are functions of rj, T* and A. 
Comparison with Equation (|16p shows that in the limit 77 — s- one must have 
A Y ->■ 0, ^ 2 0, A 3 -+ 1, B 1 - e" 1 / 71 *, 5i -> e" 1 / 7 " and £ 2 (1 - e~ 1/T ')l- 
The condition (|10p allows one to express the parameters A\, A 2 and ^3 as 
linear functions of Bq and B 2 [1(3] : 



Si 



I + 277 



l + ^-2r ? (A 3 -l) J B 2 + ^(A-l) 2 (A 2 +2A+3)y3 -£ 2 +(A-l)£ , (18) 



Qrj 1 



I + 277 



3 + 4(A a - 1)732 - (A - 1) 2 (A 2 + 2A + 3)v3 



(19) 



M= \T^^ l ~ 7] + 2{X ~ l) [ 1 - 2 ^( A + 1 )]^2-(A-1) 2 [(1-77(A + 1) 2 ]B }, 

(20) 



M = IT2^{ (1 " 77)2 " 6r?(A - x ) ( A + 1 - 2r ? A2 ) ^ 2 

+r?(A - 1) 2 [4 + 2A - r?(3A 2 + 2A + 1)]5 }- 
From Equations ((TJ) and (|17|) . we have 



(21) 



(22) 



More generally, the complete rdf follows from Equations (|12p . (| 1 3|) and (|17p . In 
particular, 



M r ) = ^w(r)Q(r) + ip u (r + l- A)6(r + 1 - A), 



(23) 



Mr) = V'2o(r)e(r) + ^ 21 (r + l-A)G(r + l-A)+V'22(r + 2-2A)e(r + 2-2A), (24) 
where 

3 



(25) 



-02fc(r) = -4vr 2 J^ 



i=l 



rW 2k {s i ) + W! 2k {s i )-W 2h {s i ) — 



A"( Si ) 



A'( Si ) 



(26) 
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Here, Sj are the three roots of A(s) = — I2r] + A\s + A 2 s 2 + A3 s 3 , the primes denote 
derivatives with respect to s and the functions Wik(s) and W 2 k(s) are defined as 

W w (s) = 1 - B + B lS , W 11 (s) = B + B 2 s, (27a) 
W 20 (s) = s[W w (s)}\ W 22 (s) = s[Wn(s)] 2 , (27b) 

W 21 (s) = 2sW w (s)W n (s). (27c) 

To close the proposal, we need to determine the parameters Bq and B 2 by im- 
posing two new conditions. An obvious condition is the continuity of the cavity 
function at r = A, what implies 

g(\-) = e- 1 / T 'g(\ + ). (28) 

Prom Equations (|12p and (|23p . assuming A < 2, one has 



g(i+) = ^-iKfe(0) (M) 



Thus, Equation (J28D yields 

ipio{\ - 1) = /T . — ■ (31) 

Equations (fT3j) and (fT7l) imply V'ii(O) = 2itB 2 / A3. Therefore, making use of Equa- 
tion ([25]) . one gets 

2^ Al + 2A 2St + 3A 3 sj 4 " (eVT* - 1) A 3 ' ld j 

Note that, taking into account (13 1 1) . Equations (129|) and fj30f) can be rewritten as 



= 1(1-^-)^ - < 34 > 

As an extra condition, we could enforce the continuity of the first derivative y'(r) 
at r = A 0| . However, this complicates the problem too much without any relevant 
gain in accuracy. In principle, it might be possible to impose consistency with a 
given equation of state, via either the virial route, the compressibility route, or 
the energy route. But this is not practical since no simple equation of state for SS 
fluids is at our disposal for wide values of density, temperature and range. As a 
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compromise between simplicity and accuracy, we fix the parameter Bo at its exact 
zero-density limit value, namely [3] 

B = l-e- 1/T '. (35) 

Therefore, Equation (I32p becomes a transcendental equation for B2 that needs to 
be solved numerically. 

In summary, Equations (fl8j) - (|2ip . (|32j) and (f35l) provide the coefficients A\, A2, 
A^, Bq, B\ and B2 as functions of 77, T* and A. This in turn determines the 
Laplace transform G(s) via Equations (|4|) and (|17j) . The rdf g(r) can be obtained 
by numerically inverting G(s) or, alternatively, by means of Equations (|12p and 

PD. 

As a further asset of our present formulation, we must point out that, once G(s) 
is determined, the static structure factor S(q) (where q is the wavevector ) of the 
SS fluid may be readily obtained as 



S(q) = l + p J dve- l ^[g{r)-l] 
G{s) - G{-s) 



1 - 2vrp 



(36) 

s=\q 



where i is the imaginary unit. 

Finally, as proven in the Appendix, the RFA proposal (I17p reduces to the exact 



solution of the PY equation 44|, [45J in three independent HS limits: (a) A — > 1, (b) 
T* -»• 00 and (c) T* -> 0. 

Before closing this section a further comment is in order. Our RFA approach has 
been based upon the exact rdf, g ex (r), at zero density [cf. Equations (|14p - (ll6p ]. On 
the other hand, it can be shown that, to first order in density, <7RFA( r ) differs from 
9e*( r ) hi the region 1 < r < I. More specifically 0], 

limMWi/T- (1 _ e -V*- W i) 

•q— >0 7] V / 

x^[(/-l) 2 -3(/ + l)(r-l)] 

x©(r - 1)0(/ -r), (37) 

where Ag(r) = <7rfa i r ) — 3sA. r )- This drawback could be corrected by a modified 
RFA (mRFA) of the form Q 

9mRFA(r) = gRFA(r) exp [-r]Q(r)Q(l - r)] (38) 

with 

Q(r) = — [Qo-Qi(r-l)], (39) 
r 

such that lim^o Qo = (l - e^ 1 /? 1 *) (z_i)3 and lim^o <5i = 3 (l - e" 1 / 71 *) (/ 2 -l). 
The coefficients Qq and Qi for finite r\ could be determined by imposing two extra 
conditions, such as the continuity of y'{r) and y"(r) at r = I. Nevertheless, since 
such conditions again lead to technical complications and their usefulness is not 
known a priori, for the sake of simplicity we restrict ourselves here to the unmodified 
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Figure 1. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.5, 
T* = 2 and rj = 0.3142 (per 3 = 0.6) as obtained from the RFA approach (solid line), the PY equation 
(dashed line), the HNC equation (dotted line) and simulation data from [40l | (circles). 

RFA described by Equations ©, ([r7|) - ([2T]) . (|52j) and (J35J). The limitations imposed 
by this restriction may be assessed from the examination of the final results. 

3. Comparison with simulation data and integral equation theories 

In order to assess the value of our theoretical approximation for the structural 
properties of SS fluids, in this section we perform a comparison between the results 
obtained with our approach and those obtained both from simulation and from our 
own numerical solution of the PY and HNC integral equations. 

As far as we know, the available simulation data for the rdf of the SS fluid as 
a function of distance at various packing fractions are those of Lang et al. 0], 
Zhou and Solana and Guillen-Escamilla et al. [4(|. Although we have made a 
comparison with all these data, in Figs. [THE] we only show graphs of g(r) vs r for 
some representative cases. Since our development is inspired by and reduces to the 
form of the solution of the PY equation for HS fluids (see the Appendix), we have 
also included in the figures the results that we have obtained from the numerical 
solution of the PY equation for the SS fluid, as well as those stemming out of the 
HNC equation for the same system 

We start by fixing a shoulder width parameter I = 1.5 and a relatively high 
reduced temperature T* = 2. Figures [T] and [2] show g(r) for packing fractions 
T] = 0.3142 and r\ = 0.4189, respectively. We observe a general good behaviour 
of the RFA, except that it underestimates the rdf near r = 1. Interestingly, the 
performance of the RFA is clearly better than that of the numerical solution of 
the PY integral equation. On the other hand, the numerical solution of the HNC 
integral equation shows an excellent agreement with the simulation data. 

The above features become much more apparent at a relatively low reduced 
temperature T* = 0.5, as Figs. [3] and d] show. At rj = 0.2094 the limitations of 



lr rhe numerical solutions were obtained by solving the system of algebraic equations resulting from the 
discretisation of the integral equation. The convergence of the solutions was found to be acceptable for a 
grid size of Ar = 0.0125 and a cut-off distance r max = 4. 
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Figure 2. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.5, 
T* = 2 and r\ = 0.4189 (po -3 = 0.8) as obtained from the RFA approach (solid line), the PY equation 
(dashed line), the HNC equation (dotted line) and simulation data from [40l | (circles). 
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Figure 3. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.5, 
T* = 0.5 and rj = 0.2094 {per 3 = 0.4) as obtained from the RFA approach (solid line), the PY equation 
(dashed line), the HNC equation (dotted line) and simulation data from [39l | (circles). 



the RFA are essentially restricted to the contact region r > 1. However, the RFA 
becomes much less reliable at the larger density rj = 0.4. Again, the PY is rather 
poor, even at rj = 0.2094, while the HNC keeps being very good. 

It is interesting to analyse the influence of temperature for a fixed value of I and 
rj. This can be achieved by comparing Figs. [2j [4] and [5] (even though the packing 
fraction in the case of Fig. [2] is not exactly the same as that of Figs. |4]and[5]). As 
the temperature decreases, the contact value increases moderately and also g(r) 
for r > 2 becomes more structured, especially when going from T* = 1 (Fig. [5]) 
to T* = 0.5 (Fig. 3]). The strongest influence of temperature occurs in the region 
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Figure 4. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.5, 
T* = 0.5 and r\ = 0.4 as obtained from the RFA approach (solid line), the PY equation (dashed line), the 
HNC equation (dotted line) and simulation data from Q (circles). 



5r 

4 \ 




1 

1.0 1.5 2.0 2.5 3.0 

r 

Figure 5. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.5, 
T* = 1 and rj = 0.4 as obtained from the RFA approach (solid line), the PY equation (dashed line), the 
HNC equation (dotted line) and simulation data from Q (circles). 



r ~ I, the discontinuity at r = I being much more pronounced as the temperature 
decreases, as expected on physical grounds. 

To assess the influence of the shoulder width we consider the choice I = 1.2 in 
Figs. [6] and [7j both at a common packing fraction r] = 0.4. As before, the RFA 
underestimates g(r) near r = 1, this effect being more important in the case of 
the PY theory. However, the HNC solution is very reliable in that region. Upon 
comparison between Figs. [5] and [6l on the one hand, and Figs. [4] and [71 on the 
other hand, one may observe that shrinking the shoulder at fixed temperature and 
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density makes the RFA and, to a lesser extent, the PY approximation become more 
reliable, while the HNC approximation is slightly less accurate. 

The above conclusion is consistent with the fact that the SS model becomes 
closer and closer to the HS model as the shoulder width decreases. In this HS limit 
the RFA reduces to Wertheim-Thiele's [4Jj exact solution of the PY equation (see 
the Appendix) and it is well known that such a solution is much more accurate 
than the HNC numerical solution for the HS fluid. Since the HS potential is also 
reached from the the SS one in the high-temperature limit (T* —¥ oo), a better 
performance of both RFA and PY over HNC can be expected to hold for sufficiently 



March 21, 2011 0:44 Molecular Physics man'subm 



12 



S. B. Yuste, A. Santos and M. Lopez de Haro 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 



Figure 8. Radial distribution function g(r) as a function of distance r for an SS fluid having A = 1.3, 
T* = 5 and r\ = 0.4189 (per 3 = 0.8) as obtained from the RFA approach (solid line), the PY equation 
(dashed line), the HNC equation (dotted line) and simulation data from [40ll (circles). 



high temperatures. This is confirmed by Fig. [8] in the case I = 1.3, T* = 5 and 
rj = 0.4189 (per 3 = 0.8). In this state, the HNC theory clearly overestimates the 
contact value, a characteristic feature of the HS system. The SS model also becomes 
equivalent to the HS model (this time with a hard-core diameter la) in the zero 
temperature limit (T* — > 0). Therefore, again the RFA and PY predictions should 
be superior to the HNC ones for sufficiently low temperatures. From a practical 
point, however, this requires temperatures so low that g(r) ~ for 1 < r < I, what 
is obviously not the case of T* = 0.5, as Figs. EJ [Hand [7] show. 

In summary, from the analysis of the results it is clear that the RFA approach 
represents a clear improvement as compared to the PY approximation. It also does 
a fair job in comparison with the HNC equation, which is generally the best ap- 
proximation. On the other hand, this latter becomes more inaccurate, particularly 
at the contact value, as the HS limit is approached, in contrast to the good perfor- 
mance of the RFA approach in this limit. Also worth mentioning is the fact that, 
although not shown, the RFA also beats the numerical solution of the OZ relation 
closed by the hard-core condition g(r) = for r < 1 plus the parametrization of the 



direct correlation function for r > 1 proposed by Guillen-Escamilla et al. 40j]. As 
expected, although the RFA always underestimates the contact value, the present 
approximation works particularly well at all distances for narrow shoulders at high 
temperatures and low densities. As the shoulder width increases, for a fixed (rel- 
atively high) temperature, the packing fraction at which deviations become more 
pronounced decreases. In the case of low temperatures, the trend observed for low 
densities is maintained but this time the region between contact and the shoulder 
width is described more poorly, while for greater distances the approach seems to 
capture rather well the subsequent oscillations of the rdf. 



4. Concluding remarks 



In this paper we have presented an (almost completely) analytical procedure to 
obtain the structural properties of SS fluids. Although the derivation heavily relied 
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upon a parallel development for the SW fluid [H, the results could not be antici- 
pated since the SW and SS potentials are physically quite different: while the SW 
potential has an attractive part that allows for the existence of a vapour-liquid 
phase transition ending at a critical point, the SS potential is purely repulsive. 

Our procedure, which follows the same rationale that we have used for other 
systems is inspired on the analytical solution of the PY equation for HS 



fluids [H, 45] (to which it reduces in the appropriate limits) and represents a good 
alternative to the usual integral equation approach of liquid state theory, which 
requires numerical work. As a matter of fact, as shown above, it clearly provides 
an improvement over the results of the simplest such integral equation for SS fluids, 
namely the OZ equation with the PY closure, and compares reasonably well with 
the results of the HNC integral equation, which in particular it beats near the HS 
limits. 

The importance of analytical or semi-analytical approximations for the equi- 
librium structural properties of simple fluids cannot be overemphasized. In this 
respect, we find it especially remarkable the fact that our approach, which only 
requires the solution of a single transcendental equation when the SS potential 
differs from the HS one (i.e., if A 7^ 1 and T* is finite), behaves better than the PY 
integral equation, whose solution in the same situation, however, involves much 
more numerical work. 

Through the comparison with the available simulation data 0, EHJ |4(| , we have 
been able to identify roughly the region in the (77, T*,l) space that these data span 
where the approximation for g(r) yields a reasonably good performance. The RFA 
approach for SS fluids leads to rather accurate results at any fluid density if the 
shoulder is sufficiently narrow (say I < 1.2), as well as for any width if the density 
is small enough (77 < 0.4). However, as the width and/or the density increase, the 
RFA predictions worsen, especially at low temperatures and in the region between 
contact and I. In any event, it is fair to insist on the advantage of the RFA quasi- 
analyticity and its relative good numerical results at low and high fluid densities 
for the SS fluid. 

Apart from the rdf, our proposal allows the immediate computation of the static 
structure factor, as shown by Equation (J36]> . We are not aware of the existence of 
simulation data for this structural property in the case of the SS interaction so a 
comparison is not feasible at this stage. We hope that our work may serve as a 
motivation to carry out such simulations. Also worth pointing out is the fact that 
the availability of the rdf allows us to obtain an approximate direct correlation 
function c(r) through inversion of the OZ equation in Fourier space. Once more, 
we are not aware of any simulation results for the c(r) of this system one could 
compare with. Given this situation and since these structural properties do not 
exhibit any particular feature we have not presented plots of them. However they 
may be easily produced upon request. 

The present results suggest that the study of the structural properties of fluids 
whose particles interact via discrete potentials composed of combinations of square 
wells and square shoulders, and other piece-wise potentials, may be tackled in a 
similar way. We plan to carry out some of these studies in future work. Finally, 
we also plan to examine in the near future the prediction and discussion of the 
thermodynamics of SS fluids which follows from our approach. 
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Appendix A. The hard-sphere limit 



In this appendix we will show that the proposal we have introduced for the structure 
of the SS fluid reduces, in the appropriate limits, to the PY approximation for HS 
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fluids. We begin with the latter. For an HS fluid, the Laplace transform G (s;rj) 
of r<7 HS (r) in the PY approximation, where we have made it explicit that it depends 
on both s and the packing fraction 77, may be expressed as in Equation (j4|), where 
the auxiliary function ^ hs (s;t/) takes a form similar to that of Equation (I17|) . 



namely |43l. 144 1 



HS 1 -127? + A* s (ri)s + A% s ( V )s 2 + A$ s (ri)s 3 

Z7r l + B l (7])s 



with 



4%) = ^, (A3) 



^»(,) = «2^, (A4) 



and 



For the sake of clarity, it is convenient in the case of the SS fluid to include 
explicitly the dependence on the packing fraction 77, the interaction range I, and the 
temperature parameter T* in the function G(s) through the notation G(s; rj, I, T*). 



A.l. Limit Z — ^ 1 

We must clearly recover the HS case if I = 1, in which case the SS potential becomes 
equivalent to an HS interaction of diameter a = 1, i.e. 

G( S ;r ? ,l,r*) = G HS ( S ;r / ). (A6) 



It is straightforward to see from Equations (119p - (121j) that in the limit I — > 1 
one has Ai(rj, l,T*) — > ^4f s (r/), irrespective of the values of lim/^i B${r), I, T*) and 
lim/^.i B2(ij, l,T*). On the other hand, the denominator of Equation (117p becomes 

TTQ 

1 + [51(77, 1,T*) + ,62(77, 1, r*)] s = 1 + B x (r])s on account of Equation JTSJ. This 



completes the proof that Equation (|17p reduces to Equation (lAlh in the limit Z — s> 1 
for arbitrary T*. 



A.2. Xirmf T* ->■ 00 

The high-temperature limit T* — ?• 00 can be understood as the limit e — >■ at finite 
T, so the SS potential ([2]) trivially becomes that of HS. Consequently, 



G(s;r,,l,oo) =G Hb (s; V ). 



(A7) 



March 21, 2011 0:44 Molecular Physics man'subm 



16 REFERENCES 

In the limit T* — > oo one has e l l T — > 1, so that Equations (I32p and ([35]) im- 
ply B Q (ri, I, oo) = B 2 (r], I, oc) = 0. Next, Equations (fTH]) -(|2T1) yield B^rj, I, oo) = 

£>! (rf) and ^(77, 2, 00) = AY s (r]) for any Z. 
A.3. Limit T* -> 

On the other hand, the low-temperature limit T* — > corresponds to the limit 
e — > 00 at finite T. In that case, the SS potential becomes equivalent to an HS 
interaction of diameter la = I. This latter condition implies that 

G(s; 7], 1, 0) = l 2 G RS (ls; I 3 ?]), (A8) 

which is a non-trivial scaling relation. In turn, taking into account the relationship 
(@D, Equation (jA8|) leads to 

*(s; V ,l,0) = ^^/ HS (/ S ;Z 3 r / ). (A9) 

From Equations (|17p and (|Aip it follows that Equation (|A9j) is satisfied if the 
coefficients comply with the following conditions: 

B ( V J,0) = 1, (AW) 

Bi(t/,I,0)=0, (All) 

Bafo.I.O) =/^f S (/ 3 r ? ), (A12) 

AifoZ.O) = l~ 2 Af s (l 3 r]), (A13) 

i 2 ( J] ,l I 0) = r 1 A« s (^), (A14) 

A 3 (r ? ,/,0)=^S(/ 3 7 ? ). (A15) 

Equation (|A10|) is a direct consequence of Equation (|35|) , Next, Equation (|32|) in 
the limit T* ->• yields Equation l|ATTj) . Then, by imposing £?i = in Equation 
(118|) one gets Equation (|A12|) . Finally, by inserting Equations (lAlOj) and (IA12|) into 
Equations () 19f ) — (12 1 j> it is easy to check that Equations (lA13j) - (|A15|) are verified. 



